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Abstract. A simple model for motion of dust particle (meteoroid) under the action of 
(solar) electromagnetic radiation is presented. The particle of the form of plane mirror 
is taken into account and exact analytical results are presented. As for long-term orbital 
evolution, particle may spiral outwards the central body (Sun); initial conditions are 
important. As a consequence, motion of real dust particles may differ from that generally 
considered. 



1. Introduction 

The Poynting-Robertson effect (P-R effect) (Robertson 1937; Klacka 1992a - the most 
complete form of the P-R effect) is generally considered to be the real effect which causes 
inspiralling of interplanetary dust particles (IDPs), meteoroids, towards the Sun. (Other, 
more simple correct derivations may be found in Klacka's papers: 1992b, 1993a, 1993b). 
The most general case of the validity of the P-R effect requires that Eq. (120) (or, Eq. 
(122) for the moving particle) in (Klacka 1992a) holds. This may not be the case for the 
real nonspherical particle, as it was discussed in Klacka (1993c, 1993d) and applied in 
Klacka (1994a). 

General equation of motion of IDP in terms of optical properties was presented by 
Klacka and Kocifaj (1994) - the paper does not present any quantitative calculation. 

The aim of this paper is to make a detailed correct calculations for plane mirror 
particle as for its equation of motion and as for secular changes of orbital elements. It 
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is considered that particle rotates around an axis fixed in space. The advantage of this 
simple model is that it can be treated in an analytical way and it may show possible 
orbital evolution of real dust particles. 

2. Equation of motion 

Let us consider a particle of a plane form, both sides of which are covered with mirror. 
The area of each side is A. When the particle's position is characterized by radius vector 
r with respect to the source of electromagnetic radiation, it's instantaneous velocity is v 
with respect to the source. 

Density of the flux of the radiation energy of the source is S when measured in the 
rest frame of the source. If c is the speed of light (all occurs in vacuum) and S = r/\r\, 
then the quantities measured in the rest frame of the particle (proper inertial frame of 
reference) are: 



where unit vector S' a represents direction and orientation of the beam reflected from the 
plane mirror characterized by normal unit vector n'; all quantities are correct to the first 
order in v/c. 

If the plane mirror is turned at an angle 5 about an axis characterized by unit vector 
a, the new normal unit vector n' n (coordinate system is fixed - it does not change) is 

n' n = n' cos S + a (a ■ n') (1 — cos 5) + (a x n') sin 5 . (2) 

We will consider that n' = ±x' = ±x and the particle rotates around axis x', for the 
sake of simplicity. The effective area for the incident radiation is 



S[ = (1 + v ■ S/c) S - v/c 
S' Q = S>-2 (n> ■ S>) n> , 
S' = (1 - 2 v ■ S/c) S , 



(1) 



A e ff = A I COS 8 



(3) 



where & is the angle between S and x (r = x x + y y + z z). We can write, thus: 



ri = - x' = - x , Q e < - tt/2, tt/2 > 



n' = +x' = + x , 6e< tt/2, 3 tt/2 > . 



(4) 



Equation of motion for the particle has the following form 




(5) 
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in the proper inertial frame of reference. 

Using the previous relations (and Lorentz transformation to the first order in v/c), wc 
finally obtain equation of motion in the reference frame of the source of electromagnetic 
radiation: 

6 = 2 — ffl - --S) S x - --x) Icoseix, (6) 
mcLV c / c J 

where m is mass of the particle and the dot denotes differentiation with respect to 
time. 

3. Electromagnetic radiation and gravitation 

Taking into account gravitational acceleration of the central body of mass M, we can 
write the complete equation of motion in the form 



^4 = — . ro 

r z m c 

where fj, = G M and G is gravitational constant. 



4. Secular changes of orbital elements 

Let the initial position and velocity vectors lie in the xy— plane. Eq. (7) yields that the 
orbital plane conserves. Secular changes of other orbital elements can be also calculated 
in an analytical way. The results are summarized in Eqs. (8) - (10) for semimajor axis a, 
eccentricity e and longitude of pericenter ui (the quantity p is defined p — a (1 — e 2 )). 

da 4 (3 fj, 2 / 3 + [6 cos (2 u) / 5 - 1] e 2 

< dt > ~^~c a (1 - e*f 2 U 



de 2 (3 JJi f 4 , , N . .„ , 



OO .tt/2 

+ (1 - e 2 ) V e 2 ™" 1 / cos 2 6 sin9 cos 2 " (9 - oj) d Q + 

n=2 J - "/* 



\JVT~p 

— e 



2 22 4 

- - — sin uj + — cos (2 w) + 

5 15 15 



2 ^ j T / 2 cos 2 6 sine cos(e — lj) sin (e — lj) 
+ - 1 ' 6 > I ff/2 1 - e 2 cos 2 (6 - d 6 



(9) 



(iw 2 /3 JJi f 2 

< 7T > = -n^~e 3 C ° S " 
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N 2n 



COS (2 lu) 

+ 



2 n + 1 2 n + 3 



+ 



+ sin(2w) sin w ^ (e sinw)' 

n=0 

p 2n /-tt/2 - w 

cos f" 

rr, ^ " + ■> ./ , , 

N 2n 



00 2n /-ir/2 - w 

(2 w) V / cos 2n+1 x d x - 

' £j 2 „ + 3 J_ „ /2 _ „ 

3 „• ^ , a ,■ 3 /o .i ( e smw ) 



— sin (2 w) + 4 sin 3 w cos (2 lu) ^ 

n— 



(10) 



Initial values for osculating orbital elements must be calculated for the case when central 
acceleration is given by the value \i - as for comparison see Klacka (1992c, 1999). 

5. Discussion to secular changes of orbital elements 

Eq. (10) shows that advance of pericenter (perihelion, in the case of Solar System) exists. 
If the relativistic terms proportional to Jifp / c would not exist, then < duj/dt > would 
be zero for cos a; = 0, i. c. for lu = ±ir/2. The relativistic terms yield that < duj/dt > < 
for lu = +7r/2 and < dco/dt > > for co = —ir/2. 

The stable situation occurs at to = —Tt/2 + e, where e is small positive number (de- 
pends on the value of eccentricity). Eq. (9) yields that < de/dt > < 0. Eq. (8) yields that 



< da/dt > < for e > ^10/33 (approximately) and < da/dt > > for e < ^10/33 



(approximately). Thus, if the initial eccentricity is greater than e > ^/10/33 (approxi- 
mately), the particle spiral towards the central body, its eccentricity still decreases. When 



e < ^/10/33 (approximately) occurs, the particle spiral outwards the central body. 
6. Conclusion 

We have shown, in an analytical way, that special type of particle exhibits orbital evo- 
lution which is not consistent with that for the Poynting-Robertson effect (moreover, 
the behaviour depends on initial conditions) . Therefore, the P-R effect cannot be simply 
applied to the study of motion of real particles. (Numerical calculations for real particle 
are in preparation.) 
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